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High  specific  impulse,  high  specific  power  devices,  such  as 
magnetoplasmadynamic  arc jets,  laser  or  microwave  propulsion 
channels,  and  MHD  generators,  involve  the  flow  of  modest 
temperature  (0.5  -  5  eV)  partially  ionized  gases  at  speeds  of 
5-20  km/sec.  The  interactions  of  such  flows  with  solid 
surfaces  containing,  channeling  or  penetrating  the  flow  provide  a 
principal  source  of  concern  for  the  efficiency  and  lifetime  of 
high  specific  power  systems.  The  present  report  describes  a 
basic  research  effort  that  examines  the  plasma/surface  inter¬ 
action  experimentally  in  an  arrangement  providing  diagnostic 
access  usually  unavailable  in  mission-oriented,  device-develop¬ 
ment  projects.  Spectroscopic  techniques  are  used  to  measure  the 
evolution  and  penetration  of  surface  ablatants  into  a  partially 
ionized  flow  from  an  argon  plasmajet.  (This  report  deals  specifi¬ 
cally  with  flow  over  a  Teflon  wedge) . 

A  theoretical  formulation  is  developed  that  allows  con¬ 
sideration  flow  in  a  regime  for  which  collisional  mean  free  paths 
(for  momentum  transfer)  are  comparable  to  scale  sizes  in  the 
flow.  A  transition  flow  shock-layer  model  is  constructed  and 
applied  to  the  wedge  flow  (in  the  manner  of  oblique  shock 
analysis) .  The  model  is  then  extended  to  include  leading  edge 
interaction  with  the  shock-layer  structure.  The  calculated  flow 
field  is  used  as  the  basis  for  estimating  penetration  of  the  flow 
by  surface  ablatants.  Comparison  of  theory  and  experiment  is 
performed  indicating  that  the  measured  ablatant  densities  in  the 
flow  are  due  to  diffusion  within  a  few  mean  free  paths  of  the 
leading  edge  (a  few  tenths  of  a  millimeter  in  the  experiments) . 
This  suggests  that  the  lifetime  of  components  in  high  specific 
power  flows  can  depend  on  the  initial  details  of  transitions  in 
flow  direction.  A  benchmark  experimental  test  stand  to  develop  a 
flow  data  base  in  the  regime  appropriate  to  the  actual  device  of 
interest  is  recommended  as  a  useful  part  of  a  lifetest  and 
evaluation  process. 
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II.  INTRODUCTION  J 


The  interaction  of  a  high  speed  (5-20  Jan/sec)  plasma  flow 
of  modest  temperature  (0.5-5  eV)  with  a  solid  surface  is  a  basic 
phenomenon  in  a  variety  of  high  specific  power  devices,  such  as 
advanced  high  specific  impulse  thrusters.  Study  of  the  details 
of  processes  involved  in  the  immediate  vicinity  of  the  surface 
is  normally  precluded  by  the  very  limited  diagnostic  access 
afforded  in  mission-oriented  devices.  The  present  research 
program  establishes  a  plasma  flow  by  means  of  a  quasi-steady 
magnetoplasmadynamic  arc jet  and  exposes  simple  solid  surfaces 
to  this  flow  while  examining  the  plasma/surface  interaction 
spectroscopically.  Detailed  measurements  provide  benchmarks 
for  theoretical  modeling  that  may  then  be  applied  to  the  more 
complex  geometries  of  actual  plasmadynamic  devices. 

The  principal  difficulty  even  in  simplified  geometry  is  that 
the  flow  regime  is  neither  continuum  nor  free  molecular.  With 
particle  densities  of  n  ■  1015  -  1016  cm”3,  and  collision  cross- 
sections  on  the  order  of  Q  ■  20  -  40  x  10”16  cm2,  mean  free  paths 
(i  -  1/nQ)  for  processes  of  interest  are  in  the  range  of  0.1  to 
5  mm.  The  viscosity  of  the  flow  is  on  the  order  of  2  x  10“4 
kg/m-s  (e.g.,  argon  at  -  1  eV)  with  a  mass  density  of  about 
3  x  10”4  kg/m3  and  a  speed  of  104  m/s,  so  the  characteristic 
dimension  for  viscous  effects  (defined  by  Rey  ■  p -  1)  is 
L  »  0.07  mm  «  i.  Thus,  interactions  of  the  flow  with  solid 
boundaries  are  in  the  so-called  transition  flow  regime. 
Examination  of  these  interactions  therefore  must  proceed 
theoretically,  without  the  benefit  of  well-established  theories 
(particularly  if  multiple  species  and  states  are  included) ,  and 
experimentally  on  a  scale-size  that  requires  inspection  at  sub¬ 
millimeter  distances  from  the  solid  surface.  As  a  further 
complication,  both  theoretically  and  experimentally,  the  outer 
flow  field  near  the  surface  sample  is  not  easily  characterized 
since  it  involves  a  finite  size  object  in  a  highly  supersonic 
flow.  For  application  of  the  plasma  flow  normal  to  a  flat 
surface,  the  flow  resembles  that  around  a  blunt  body.  Use  of  a 
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wedge  (prism  shape)  provides  a  more  tractable  flow  field  for 
theoretical  modeling  but  still  can  suffer  from  edge  effects 
•xperimentally.  To  attack  the  plasma/ surf ace  interaction  in  the 
face  of  these  difficulties,  we  proceed  with  detailed  spectro¬ 
scopic  examination  of  the  distribution  of  species  near  the  sample 
surface,  followed  by  construction  of  simple  kinetic  models  to 
understand  the  experimental  results  and  offer  guidance  for 
interactions  in  more  complex  geometries. 
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III.  EXPERIMENTAL  APPARATUS 

The  present  experiments  have  been  performed  with  facilities 
and  diagnostics  improved  from  previous  plasma/ surface  interaction 
experiments.  The  power  source  was  upgraded  from  a  20  kJ,  150 
n sec  pulse  forming  network  to  a  400  kJ,  0.5-1  msec  LC-ladder 
style  pulse  line.  Appropriate  changes  in  the  current  transmis¬ 
sion  lines  and  coupling  to  the  arc  jet  chamber  and  in  the  arc  jet 
electrical  feed  were  also  accomplished.  The  arc jet  feed  design 
allows  several  types  of  arc jets  to  be  tested  without  changing  the 
basic  power  and  gas  handling  connections.  In  addition,  the  gas 
feed  was  upgraded  to  match  better  the  higher  currents  (75  kA  vs 
12  kA)  provided  by  the  new  power  source.  Figure  1  shows  the 
arrangement  and  components  of  the  arc jet  facility  at  RDA/WRL. 
Voltage,  current  and  electrical  probe  data  acquisition  has  been 
upgraded  to  a  computer  controlled  digital  data  acquisition  in  the 
screen  room.  The  optical  diagnostics  has  also  been  upgraded  to  a 
computer  controlled  digital  spectroscopic  data  acquisition 
system  that  allow  simultaneous  observation  spectrally  and 
spatially.  System  upgrades  have  had  the  benefit  of  activities  in 
closely-related  programs  at  RDA/WRL  on  dense  plasma  jet  propaga¬ 
tion  for  SDIO/IST  (managed  by  AFOSR/NP)  and  MPD  arc jet  thruster 
studies  for  AFOSR/NA,  including  a  cooperative  effort  with  the 
AFOSR-sponsored  MIT  group  under  Prof.  M.  Martinez -Sanchez. 

An  existing  400  kJ  capacitor  bank  was  modified  into  an  eight 
stage  LC-ladder  network  to  provide  a  flat  quasi-steady  current 
pulse  to  the  arejet.  Each  stage  is  made  up  of  16  20  kv  capaci¬ 
tors,  15.6  pF  each,  connected  in  parallel  giving  a  total 
capacitance  of  2.0  mF.  A  schematic  of  the  eight  stage  LC-ladder 
is  shown  in  Figure  2.  Each  section  is  connected  by  a  10.1  /ih 
solenoid  inductor  with  an  approximate  parasitic  resistance  of 
about  three  milliohms.  The  final  solenoid  inductor  is  20  *ih  and 
the  bank  matching  resistor  is  120  milliohms.  The  pulse  forming 
network  is  designed  with  a  high  voltage  ignitron  switch  at  each 
stage  so  that  the  output  waveform  can  be  modified  depending  on 
the  sequencing  of  ignitrons.  If  ignitrons  #1  and  #8  are  fired 
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simultaneously,  yielding  the  circuit  shown  schematically  in 
Figure  2,  the  current  pulse  duration  is  approximately  500  Msec 
with  a  rise  time  of  100  hubc.  The  pulse  forming  network  is 
connected  to  the  arc jet  by  ten  RG  213/u  coaxial  cables,  sixty 
feet  long,  with  the  center  conductor  hot  and  negative. 

The  arc jet  assembly  with  power  and  gas  feeds  is  mounted  in  a 
cylindrical  stainless  steel  vacuum  tank  6  meters  long  and  0.6 
meters  in  diameter.  The  vacuum  is  maintained  at  a  static  level 
of  less  than  10"5  torr  by  two  six-inch  diffusion  pumps.  The  tank 
and  associated  pumps  are  electrically  floating  to  prevent 
alternative  current  paths  for  the  PFN-arcjet  circuit.  The  design 
of  the  arc jet  plasma  source  for  the  surface  plasma  interaction 
study  is  based  on  experience  with  similar  systems  over  the  last 
twenty  years.  In  particular,  it  was  recognized  that  rather 
simple  arrangements  would  provide  adequate  performance  without 
delaying  experimental  progress  to  deal  with  long  term  concerns 
such  as  insulator  or  electrode  erosion.  Accordingly,  the 
electrodes  were  made  of  brass  and  the  insulator  (backplate) 
through  which  gas  is  injected  was  made  of  Plexiglas.  Figure  3 
provides  a  schematic  view  of  the  arc jet.  The  electrodes  assemble 
into  a  nested  conical  arrangement  that  allows  the  diffuse  arc 
processes  (largely  in  the  snubber  region  off  the  backplate)  to 
accelerate  the  injected  gas  and  provide  a  spatially  uniform  exit 
flow.  Gas  J.s  provided  to  the  arcjet  from  a  large  gas  reservoir 
placed  close  to  the  insulator  backplate.  Six  valves  feed  a  pre¬ 
expansion  chamber  directly  upstream  of  12  choked  orifices  in  the 
insulator  backplate. 

All  data  for  the  present  experiment  was  collected  with  a 
consistent  set  of  operating  conditions.  The  pulse  forming 
network  was  charged  to  13  kV,  the  initial  system  vacuum  was 
2  x  10”5  torr  or  less,  and  the  gas  reservoir  was  filled  with  300 
torr  absolute  of  argon-1.5  percent  hydrogen  mixture.  This  gas 
reservoir  pressure  yields  a  mass  flow  rate  of  6  g/s  assuming 
choked  flow  at  the  insulator  backplate  orifices.  For  all 
experiments  the  valve  drive  system  was  fired  10  ms  before  firing 
the  bank.  A  block  diagram  of  the  arcjet  system  is  shown  in 
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Figure  4.  The  arc jet  breech  voltage  was  measured  by  a  1000:1 
Tektronix  voltage  probe  optically  isolated  from  the  data 
acquisition  system  at  the  power  feed  external  to  the  vacuum  tank. 
The  total  current  was  measured  by  a  large  Rogowski  loop  surround¬ 
ing  the  cathode  sleeve  at  the  base  of  the  power  connection.  This 
probe  measures  the  time  rate  of  change  of  the  enclosed  magnetic 
flux,  which  is  integrated  to  give  the  enclosed  current.  An 
active  integrator  with  a  377  psec  integration  time  constant  and  a 
70  msec  associated  droop  time  was  used  with  the  Rogowski  loop. 

The  total  current  and  breech  voltage  were  monitored  every 
discharge  by  the  computer  controlled  data  acquisition  system. 
Figure  5  shows  the  total  current,  breech  voltage  and  delivered 
power  for  a  typical  discharge. 

III. A.  SPECTROSCOPIC  ANALYSIS  SYSTEM 

To  enable  spectroscopic  analysis  of  the  entire  flow  field, 
a  digital  spectral  and  spatial  data  acquisition  system  has  been 
developed  in  a  cooperative  effort  with  the  AFOSR- sponsored  MIT 
group  under  Prof.  Martinez-Sanchez .  The  system  is  comprised  of 
a  1.2  meter,  f/11.5  spectrograph  with  an  EG&G/PARC  intensified 
optical  multichannel  analyzer  (SIT  camera) .  The  spectral  and 
spatial  output  from  the  camera  is  processed  by  an  IBM  System  9000 
laboratory  computer  for  immediate  background  subtraction  and 
signal  averaging.  The  digital  spectral  and  spatial  data 
acquisition  system  development  consists  of  interfacing  systems, 
characterization  of  SIT  camera  performance,  relative  spectral 
calibration,  optical  alignment  with  spatial  calibration, 
development  of  gate  operation  capability  and  computer  analysis  of 
data  with  associated  software  development.  Also,  considerable 
effort  was  necessary  to  shield  the  SIT  camera  and  associated 
electronics  from  electromagnetic  noise  generated  by  the  arc jet 
power  supply  and  connections.  Many  issues  of  SIT  system 
performance  surfaced  during  system  characterization  and  calibra¬ 
tion.  These  include  linearity  of  signal,  spatial  uniformity  of 
responsivity,  background  noise  due  to  internal  electronics,  shot 
noise  due  to  spurious  electrons  from  the  photocathode,  and  edge 
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Block  diagram  of  the  MPD  arc jet  system.  ; 
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of  field  effects  that  cause  'pincushion',  distortion  at  frame 
edges  (straight  lines  are  observed  as  curves) . 

The  SIT  camera  has  a  detector  array  of  512  x  512  pixels 
(262k  pixels)  with  spatial  size  of  25  microns  and  a  dynamic 
range  of  16k.  The  interfaced  system  had  resolution  less  than 
expected  due  to  signal  blooming  from  pixel  to  pixel  (thus 
decreasing  the  observed  spatial  resolution)  and  range  less  than 
expected  due  to  camera  electronics  (e.g.,  A/D  converters)  and  to 
software  capabilities.  The  detector  system  was  therefore  limited 
to  250  horizontal  pixels  (spectral  data)  and  to  128  vertical 
pixels  (spatial  data)  or  32k  data  points  per  frame.  The  detector 
showed  excellent  spatial  uniformity  of  responsivity  but  limited 
range  of  signal  linearity  of  0  to  600  counts  as  shown  in  Figure 
6.  The  detector  system  relative  spectral  response  was  measured 
with  a  calibrated  tungsten  filament  lamp.  The  calibration  source 
was  an  Optronic  Laboratories  Model  550C  tungsten- filament  lamp. 
The  lamp  was  operated  at  15  amperes  (±  0.01  amperes)  DC  for  which 
tables  of  spectral  radiance  were  provided  and  traceable  to  NBS. 
The  uncertainty  in  current  contributed  less  than  one  percent 
error  in  the  spectral  radiance.  The  lamp  was  mounted  to  include 
effects  of  the  input  optics,  turning  mirrors  and  attenuation  due 
to  the  chamber  window  material.  The  resultant  relative  spectral 
response  normalized  to  1.0  at  700  nm  is  shown  in  Figure  7.  Two 
features  are  worth  noting.  The  system  has  a  rapidly  varying 
response  around  600  nm  and  furthermore,  below  370  nm,  the 
response  drops  dramatically  without  quartz  optics. 

III.B.  MEASUREMENTS  NEAR  SAMPLE  SURFACE 

The  experimental  task  has  been  the  detection  and  charac¬ 
terization  of  sample  species  in  the  plasma  flow  and  the  measure¬ 
ment  of  temperature  and  density  with  fine  spatial  resolution  in 
the  immediate  neighborhood  of  surfaces  exposed  to  the  arc jet 
plume.  A  schematic  of  the  typical  optical  arrangement  is  6hown 
in  Figure  8.  The  sample  to  be  investigated  is  mounted  down¬ 
stream  of  the  arc jet  exit  in  the  plume  which  affords  easy 
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Figure  6.  SIT  camera  observed  counts  (vertical  axis)  xft  relative 
input  light  flux  (horizontal  axis)  at  two  horizontal 
rows  shows  good  spatial  uniformity  of  responsivity  and 
linear  response  to  600  counts. 
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Figure  7.  Relative  system  response  yi  wavelength  normalized  to 
1.0  at  700  nm.  Quartz  optics  (* — *)  and  plexiglas 
window  (o—o)  shows  nonuniform  responsivity. 
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Figure  8.  Schematic  of  typical  optical  arrangement  for  plasma  surface 

interaction  experiments.  The  first  turning  mirror  represents  three 
first  surface  mirrors  to  rotate  desired  spatial  view  to  align  with 
spectrograph  slit. 
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optical  access  and  absence  of  electric  and  magnetic  fields  that 

would  further  complicate  the  observed  phenomena.  A  schematic 

view  of  the  Teflon  sample  used  in  the  current  experiments  is  ; 

shown  in  Figure  9.  The  sample  is  a  22  degree  wedge  with  the  a 

flat  side  on  the  arc jet  centerline  and  a  width  of  1  cm.  The 

optical  system  allowed  measurement  of  optical  emission  orthogonal  ! 

to  the  wedge  surface  with  a  spatial  resolution  of  52  microns.  j 

The  expected  observed  species  were  argon,  carbon,  fluorine  and 

hydrogen  (which  was  added  to  the  argon  mass  injection  for  j 

diagnostic  purposes) .  Spatially  resolved  spectra  of  emission  i 

from  these  species  was  performed  orthogonal  to  the  wedge  surface 

at  positions  2.5,  5,  10  mm  from  the  wedge  vertex  (see  Figure  9). 

Figure  10  shows  a  visual  representation  of  SIT  camera  data  taken 
at  the  position  2.5  mm  from  the  wedge  vertex.  The  horizontal 
axis  is  wavelength  and  represents  250  data  points  or  8.65  nm. 

The  vertical  axis  is  spatial  position  and  represents  128  data 
points  or  6.656  mm.  The  dark  band  is  the  wedge  which  apparently 
contributes  a  continuum  to  the  observed  emission.  The  large 
feature  is  Hq  at  696.28  nm  and  the  two  small  features  that  only 
appear  near  the  sample  surface  are  emission  from  carbon  ions  (C 
II)  at  658.285  nm  and  657.803  nm  (larger  feature).  Figure  11 
shows  a  intensity  scan  of  Figure  10  at  row  58  which  represents  a 
position  0.5  mm  from  the  sample  surface  and  shows  the  three 
features  described  above. 

A  complete  series  of  data  was  taken  at  the  three  wedge 
locations.  Five  argon  neutral  emission  lines  were  chosen  for 
observation  due  to  their  close  proximity  in  wavelength  thus 
reducing  relative  system  calibration  errors,  absence  of  other 
emitting  species  nearby  and  the  widest  possible  separation  of 
energy  levels  for  the  lines  observed.  The  wavelength,  energy 
level,  statistical  weight  and  decay  constant  of  the  observed 
argon  lines  is  shown  in  Table  I.  Each  spatially  resolved  argon 
emission  line  was  zero-corrected  and  integrated  spectrally  and 
corrected  for  relative  system  responsivity  to  yield  relative 
emission  versus  spatial  position  in  relation  to  the  sample 
surface.  In  addition,  3  to  6  shots  at  each  wavelength  location 


Figure  9.  PI a ana /surface  interaction  geometry  is  a  22*  Teflon  wedge  mounted 
5.0  cm  downstream  of  arcjet.  Three  wedge  positions  orthogonal  to 
the  surface  were  investigated  as  indicated. 


visual  representation  of  emission  near  sample  surface.  Raw  data  as 
obtained  by  SIT  camera  post  shot.  Horizontal  axis  is  wavelength 
and  vertical  axis  is  position.  The  two  small  features  are  C  II 
emission  near  sample  surface  (dark  band) . 


with  ZS0  points. 
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were  averaged  in  order  to  increase  signal-to-noise  and  reduce 
possibility  of  shot  to  shot  variation.  Figures  12  shows  results 
for  a  position  2.5  nun  from  the  wedge  vertex  for  four  of  the  lines 
listed  in  Table  I.  Results  for  737.212  nm  Ar  I  were  not  included 
because  the  signal-to-noise  level  at  this  position  did  not 
provide  reliable  results.  Figures  13  and  14  show  relative 
emission  versus  spatial  position  for  the  five  argon  emission 
lines  chosen  at  positions  5  mm  and  10  mm  from  the  wedge  vertex, 
respectively. 


Table  I. 


Species 

Wavelength 

nm 

Energy  Level 
(eV) 

G 

A(Xl08) 

(1/8) 

Argon 

706.722 

5 

0.395 

738.398 

13.3 

5 

0.087 

737.212 

13.324 

9 

0.02 

703.025 

14.835 

5 

0.0278 

603.215 

15.127 

9 

0.0246 

Fluorine 

738.745 

14.763 

4 

0.38 

Carbon 

657.803 

16.35 

4 

0.48 

As  displayed  in  Figures  10  and  11,  carbon  can  be  observed  in 
emission  near  the  sample  surface  as  well  as  fluorine.  Figures 
15  -  17  show  zero  corrected,  spectrally-integrated  spatially- 
resolved  relative  intensity  for  657.803  C  II  and  703.745  F  I  at 
the  three  wedge  positions.  Table  I  shows  the  energy  level, 
statistical  weight  and  decay  constant  for  the  observed  carbon 
and  fluorine  lines.  Emission  of  sample  species  was  only  observed 
near  the  sample  at  all  observed  positions.  Sample  species 
emission  was  small  at  the  surface  with  peak  emission  at  a 
position  0.5  -  2  mm  from  the  wedge  and  tending  to  no  emission  at 
larger  distances  (  >2.5  mm) . 
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Figure  12.  Relative  observed  spatial  intensity  of  argon  vs 

distance  from  sample  at  a  wedge  position  of  2.5  mm 
from  vertex  at  a)  706.722  nm,  b)  738.398  nm, 
c)  703.025  nm,  and  d)  603.215  nm. 
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Figure  12.  (continued) 
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Figure  13.  Relative  observed  spatial  intensity  of  argon  vs 
distance  from  sample  at  a  wedge  position  of  5.0 
from  vertex  at  a)  706.722  nm,  b)  738.398  nm, 
c)  737.025  nm,  d)  703.025  nm  and  e)  603.215  nm. 
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Figure  14. 


Relative  observed  spatial  intensity  of  argon  vs 
distance  from  sample  at  a  wedge  position  of  10  mm 
from  vertex  at  a)  706.722  nm,  b)  738.398  nm, 
c)  737.025  nm,  d)  703.025  nm,  and  e)  603.215  nm. 
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Figure  15. 
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Relative  observed  spatial  intensity  of  a)  carbon 
(C  II)  667.803  nm  and  b)  fluorine  (F  I)  738.749  nm 
YS  distance  from  sample  at  a  wedge  position  of 
2.5  mm  from  vertex. 
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III.C  EXPERIMENTAL  RESULTS 

As  discussed ,  relative  argon  neutral  emission  was  observed 
at  the  three  sample  positions.  The  observed  intensity  is  related 
to  the  excited  state  density  by 

I*  -  n*A  (hc/A)  G 

where 

n*  -  excited  state  number  density 

A  -  radiative  decay  probability  or  Einstein  coefficient 
a  -  emission  wavelength 

G  -  geometry  and  system  sensitivity  factors 

If  the  plasma  is  in  local  thermodynamic  equilibrium  then  the 
two  most  important  quantities  that  determine  the  excited  state 
densities,  n* ,  is  the  free  electron  density  and  the  electron 
temperature.  In  LTE,  the  excited  state  density  is  related  to  the 
electron  temperature  by 

_  _*  -E/kTe 
*  npg  e 

"  "  Up(Te) 

where 

np  -  number  density  of  atoms  in  ionization  stage  p 
up(Te)  •  partition  function  or 

Op(T.)  -  |9i  .”Ei/lcTe 

The  plasma  electron  temperature  can  then  be  determined  using  the 
Boltzmann  plot  technique  by  combining  the  above  equations  and 
taking  the  natural  log. 

fiiAi'l  Ei 

in  lilAlJ  +  Constant  -  - 

Only  a  relative  calibration  of  the  detection  system  is  necessary 
(Figure  9)  to  determine  the  term  in  parentheses  and  Te  is  found 
by  least  squares  straight  line  fit  through  the  observed  data 
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points.  Figure  18  shows  a  typical  rssult  for  the  argon  lines  of 
Table  I  at  wedge  position  1  cm  from  the  vertex  and  2.6  mm  from 
the  surface.  A  least  squares  fit  to  the  data  yields  a  measured 
temperature  of  1.1  eV  with  an  error  of  +  0.18  eV  and  -  0.135  eV 
at  this  location.  The  plasma  electron  temperature  was  calculated 
by  this  method  at  25  locations  0.102  mm  apart  for  each  of  the 
three  wedge  positions.  These  results  are  shown  in  Figures  19,  20 
and  21  for  the  three  wedge  positions  2.5,  5  and  10  mm,  respec¬ 
tively.  Generally,  the  error  was  less  in  regions  where  the 
electron  temperature  was  not  changing  rapidly,  i.e.,  in  the  free 
stream,  and  was  on  the  order  of  +  20%  and  -  15%.  Near  the  wedge 
surface,  and  around  regions  of  electron  temperature  fluctuations 
the  error  was  larger  on  the  order  of  +  35%  and  -  22%  implying 
tendency  to  non-LTE  behavior  in  these  regions.  Unfortunately, 
implications  of  non-LTE  behavior  are  observed  generally  in  all 
regions  near  the  sample. 

The  local  electron  density  can  be  determined  by  measuring 
the  Stark  broadening  profiles  of  the  observed  emission  lines. 
Spatially  resolved  line  profile  measurements  were  made  of  the 
603.215  Ar  I  line  and  compared  with  line  profile  measurements  of 
Hq.  Figure  22  shows  an  axonometric  plot  of  the  spatially  and 
spectrally  resolved  Hq  emission  at  the  wedge  position  of  2.5  mm. 
Figure  23  shows  a  equal-intensity  contour  plot  of  the  same  line. 
Slight  variation  in  the  observed  profile  is  noted  at  rows  6  to  16 
(300  -  830  microns) .  If  the  ion  temperature  is  on  the  order  of 
or  larger  than  the  electron  temperature  then  line  broadening  due 
to  the  ion  temperature  (Doppler  broadening)  may  be  significant 
compared  to  Stark  broadening.  The  Doppler  broadened  profile  is  a 
Gaussian  and  the  full  width  at  half  maximum  (FWHM)  is  given  by 

AA0  ■  7.717  X  10~4  A 0  yfiTtT 

where 

A0  -  wavelength  at  line  center  (nm) 

T^  -  emitter  temperature  (eV) 

n  -  atomic  weight  in  atomic  units 

For  Ti  *  Te  «  1  eV,  the  Doppler  width  for  hydrogen  is  0.0506  nm 


Temperature  (eV)  2  Temperature  (eV 


i  20.  Calculated  electron  temperature  from  Boltzmann 

plots  YS.  distance  from  sample  at  a  wedge  position 
of  5.0  mm  from  the  vertex. 


Figure  21.  Calculated  electron  temperature  from  Boltzmann 

plots  v£  distance  from  sample  at  a  wedge  position 
of  10  mm  from  the  vertex. 
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whereas  for  argon  it  is  0.00736  nm.  Doppler  broadening  for  argon 
is  insignificant  (relative  to  observed  line  widths)  for  ion 
temperature  less  than  approximately  40  eV.  The  simultaneous 
presence  of  Doppler  and  Stark  broadening  results  in  a  profile 
known  as  the  Voight  profile.  The  convolution  is  given  by 

V(AA*)  ■  J Dop ( A  A  * )  Sta(AA  -  AA*)  d(AA*) 

where  Dop()  «  Doppler  profile  and  Sta()  *  Stark  profile.  One 
analytic  result  that  can  be  noted  is  that  if  either  broadening 
mechanism  independently  has  a  FWHM  less  than  one-third  of  the 
other,  the  observed  line  is  only  affected  by  15-20%. 

Figure  24  shows  a  Lorentzian  profile  fit  to  the  observed  HQ 
at  a  wedge  position  of  2.5  mm  and  780  microns  from  the  surface 
with  a  FWHM  of  0.164  nm  ±  0.007  nm.  If  no  Doppler  broadening  of 
the  observed  profile  is  assumed  then  the  electron  density  is 
1.11  x  1016/cxa3*  The  corresponding  Stark  width  measurement  of 
the  603.215  nm  argon  line  yields  a  electron  density  of  0.62  x 
1016/cm3.  If  the  difference  in  the  observed  electron  densities 
is  due  to  Doppler  broadening  of  hydrogen,  then  the  deconvolution 
of  the  observed  Hq  profile  (assuming  the  Stark  contribution  to 
the  smaller  electron  density)  yields  a  Doppler  width  correspond¬ 
ing  to  a  hydrogen  temperature  of  approximately  1  eV.  However, 
the  electron  density  calculated  from  argon  profile  measurements 
has  a  large  error  due  to  limitations  of  wavelength  resolution. 
Electron  density  calculated  from  the  Hq  profiles  is  within  this 
error.  The  observed  differences  in  calculated  electron  density 
between  Ha  profiles  and  argon  profiles  can  be  ascribed  to  error 
in  the  data  acquisition  system  and  to  errors  in  electron  tempera¬ 
ture  (which  enters  into  the  Stark  width  calculation)  without 
contributions  from  Doppler  broadening.  The  observed  electron 
density  versus  spatial  position  at  the  three  wedge  positions  are 
shown  in  Figures  25,  26  and  27.  The  results  at  the  three 
positions  are  very  similar  and  without  any  significant  features. 

Thus,  electron  temperature  and  electron  density  measurements 
at  three  locations  on  the  Teflon  wedge  indicate  a  quasi-steady 


Distance  from  sample  (mm) 


Figure  26.  Calculated  electron  density  vs  distance  from  sample 
at  a  wedge  location  5.0  mm  from  vertex. 
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Figure  27.  Calculated  electron  density  v§  distance  from  sample 
at  a  wedge  location  10  mm  from  vertex. 
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electron  density  on  the  order  of  l  x  l016/cm3  and  an  electron 
temperature  on  the  order  of  1.1  eV.  The  relative  line  Inten¬ 
sities  of  argon,  fluorine  and  carbon  show  strong  spatial 
components  Implying  changes  In  populations  spatially  at  the  three 
wedge  positions. 

In  order  to  qualitatively  evaluate  these  variations  In 
spatial  population  of  the  different  species  a  series  of  calcula¬ 
tions  was  made  assuming  that  the  ground  state  of  the  emitting 
specie  is  in  Boltzmann  equilibrium  with  the  emitting  state  and 
that  application  of  Sadia  equation  to  determine  neutral  and  ion 
densities  of  each  specie  is  valid.  Thus,  with  knowledge  of  Te 
and  electron  density,  the  local  population  of  argon  and  argon 
ions,  carbon  and  carbon  ions,  and  fluorine  and  fluorine  ions  can 
be  plotted  versus  distance  from  the  sample  for  each  of  the  three 
observed  wedge  positions.  Figures  a,  b,  and  c  of  28,  29,  and  30 
show  the  results  of  these  calulations.  Note  that  the  calculation 
of  each  specie  population  is  an  equation  of  the  form 

N  -  const  *  i  *  e&/Te 

where  the  exponential  term  is  large  and  therefore  small  changes 
in  Te  make  large  changes  in  calculated  population.  The  error  in 
calculation  is  given  by 

N  -  N  ±  N  |  ±exp  +  1  | 

where  «  is  the  error  in  the  electron  temperature.  For  example, 
in  calculation  of  argon  neutral  population  (Aro)  an  error  of 
0.1  eV  at  Te  *  1.0  eV  implies  an  error  of  +3.38  x  Aro  to 
-0.7  X  Aro- 

The  results  of  the  calculation  imply  that  the  total  argon 
density  is  relatively  flat  spatially  and  approximately  the  same 
value  at  the  three  wedge  positions.  The  situation  for  fluorine 
and  carbon  is  entirely  different  (as  expected).  At  the  2.5  mm 
position,  the  fluorine  and  carbon  density  are  about  constant  and 
extend  approximately  1  mm  from  the  surface.  The  nominal  carbon 
density  is  approximately  10  times  the  nominal  fluorine  density. 


Figure  28. 


Calculated  specie  populations  v&  distance  from  sample 
at  a  wedge  location  2.5  mm  from  vertex  where  symbols 
represent  (— ■)  total  of  specie,  (-  -)  neutrals,  and 
(+”+)  ions,  in  a)  argon,  b)  carbon,  and  c)  fluorine. 
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Also  nots  that  ths  carbon  and  fluorine  neutral  densities  are 
approximately  equal  and  represent  about  5%  of  the  species 
present.  At  5  am  from  the  wedqe  vertex  the  fluorine  spatial 
population  is  similiar  to  that  at  2.5  mm  but  spread  spatially. 

The  total  carbon  population  decreases  and  appears  to  become 
separated  from  the  surface.  At  10  mm  from  the  vertex  the 
fluorine  population  now  has  a  definite  feature  at  0.8  to  1.4  mm 
from  the  wedge.  The  carbon  population  is  (within  detection 
limits)  not  present  at  the  wedge  surface  with  a  feature  at  0.9  to 
1.3  mm.  The  densities  of  the  wedge  species  do  not  decrease  in 
simple  exponential  fashion  with  distance  from  the  surface,  but 
may  instead  depend  on  the  flow  field  and  to  plasma/ surface 
interactions  upstream  of  the  field  of  view. 
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IV.  THEORETICAL  FORMULATION 

In  the  low  density,  high  speed  flows  characteristic  of  high 
specific  impulse  and  high  specific  power  channels,  the  mean  free 
path  for  collisions  between  flow  particles  (atoms,  molecules, 
ions,  electrons)  can  be  comparable  to  dimensions  that  scale 
convective  and/or  diffusive  processes  in  the  flow  channel.  Such 
dimensions  may  be  the  viscous  boundary  layer  thickness,  radii  of 
curvature  for  changes  in  flow  direction,  and  various  scale  sizes 
for  electrical  current  distribution,  propellant  mixing,  and  heat 
transfer.  If  the  mean  free  paths  for  collisional  interchange  of 
momentum  and  energy  are  not  much  smaller  than  the  dimensional 
scale  of  interest,  then  it  is  difficult  to  specify  flow  velocity, 
pressure,  and  temperature  in  a  consistent  and  useful  fashion. 
Essentially,  more  detailed  information  about  the  distribution  of 
particle  positions  and  velocities  must  be  carried  through 
fundamental  physical  equations  than  is  tractable  with  existing 
mathematical  tools.  The  simplifications  represented  in  the 
development  of  classical  hydrodynamics  are  no  longer  available. 
For  example,  a  Maxwellian  distribution  of  velocities  about  a  mean 
velocity  characterizes  quantities  such  as  flow  speed,  tempera¬ 
ture,  and  pressure.  Fluxes  of  mass,  momentum,  and  energy  are 
then  computed  from  first  order  derivatives  of  these  quantities. 
Without  a  simple  distribution  function  allowing  meaningful  use  of 
average  quantities,  each  process  (e.g.,  heat  transfer  to  surface, 
force  on  a  plate)  must  be  calculated  separately  in  terms  of  the 
local  conditions  of  particle  velocities,  spatial  distributions 
and  geometry. 

Collisions  between  particles  can  alter  the  distribution  of 
particles  in  phase  space  such  that  moments  of  the  distribution 
function  change  too  drastically  to  be  usefully  related  (i.e.,  in 
the  manner  of  thermodynamic  state  variables  and  their  deriva¬ 
tives)  .  Continuum  flow  is  based  on  sufficient  collisionality 
("molecular  chaos")  to  maintain  an  equilibrium  distribution 
during  any  flow  change.  In  the  opposite  limit  of  very  low 
density  flows,  the  interactions  can  be  limited  to  known  distribu- 
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tlons  of  particles  In  ^ase  apace  colliding  only  with  solid 
boundaries  (or  fields) .  So-called  free  molecular  flow  Is  thus 
based  essentially  on  mapping  particle  trajectories  from  one 
region  of  phase  space  to  another.  Interactions  with  surfaces  are 
characterized  by  seml-emplrlcal  coefficients  that  attempt  to  mix 
simple  geometric  transfer  of  particle  paths  with  re-speclflcatlon 
of  some  portion  of  the  particle  distribution  function.  The 
Intermediate  or  transition  flow  regime  requires  calculation  of 
changes  in  the  particle  density  and  velocity  distributions  due 
to  a  limited  number  of  collisions  in  the  region  of  interest.  It 
resembles  free  molecular  flow  in  that  particle  distributions  are 
specified  at  least  in  some  portion  of  the  flow,  and  accommoda¬ 
tions  with  surfaces  must  be  invoiced.  Collisional  events  between 
particles,  however,  are  included  so  the  flow  will  depend  in  part 
on  atomic/  molecular  or  ionic  structure.  The  principal  diffi¬ 
culty  is  integration  of  collision  events  over  particle  veloci¬ 
ties,  especially  for  particles  that  have  already  had  the 
opportunity  for  collision.  Some  progress  might  be  possible, 
however,  by  ignoring  most  details  in  favor  of  average  particle 
velocities.  The  following  analyses,  based  on  such  average 
particles,  develop  models  in  the  transition  flow  regime  for  the 
deflection  of  a  low  density  flow  and  its  penetration  by  surface 
ablatants . 

FORMULATION  FOR  TRANSITION-FLOW  SHOCK-LAYER 

Consider  a  flow  of  atoms  with  initial  density  nQ,  as 
sketched  in  Figure  31.  Treated  as  a  beam-like  flow,  the  number 
density  n^  of  this  flow  will  be  diminished  by  collisions  with  the 
slower-moving  atoms  of  density  ns: 

d°A 

UA  dx~  *  ”nAnsQ (UA  "  us) 

where  u&,  us  are  the  mean  speeds  of  the  original  atomic  flow 
and  the  slower  flow  atoms,  respectively.  (Note  that  with 
identical  atoms  in  both  flows,  we  are  not  attempting  to  ascribe  a 
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depletion  of  the  original  flow  unless  u8  <  uA,  regardless  of  the 
contributions  to  the  collision  frequency  by  random  motions  in 
either  flow. )  Similarly,  the  growth  of  the  slow  flow  number 
density  is: 


Thus, 


dn8 

u»  djT  “  nA«sQ(UA  “  u8) 
d»A  dn8 

"a  sr  +  u»  sr  - 0 


so  nAuA  +  n8u8  ■  constant  *  n0uA. 

The  preceding  relationship  expresses  continuity  of  mass  in 
the  total  flow.  Conservation  of  momentum  and  energy  will  be 
examined  later. 


With 


n8  -  (nD  -  nA)UA/u8 
uA  >  m,  we  have: 


UA  dx~  m  “nA<no  “  nA)Q(uA  -  us)M 


eGT  -  -  nA  (nc  -  nA)  Q  (M  -  1) 
Let  nA  «  »»An0,  and  x  -  A/n0Q,  then 
d»;  A 

no2Q  37"  -  -nA(l  -  »?A)n02Q(M  -  1) 


Thus, 


“^aC1  “  »?a)  (m  “  1) 
"»?a(1  “  ha)#1  I 


»»aU-»»a) 


-ft  dA 


With  U  rn  M  -  1. 
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where  n r  and  a r  are  reference  values, 
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From  Dwight 


i,  101.1): 


_ £& _  »  .1  loa  A±feS 

x(a+bx)  a  Aog  a 


f  n  x  m  log  •  «  “  (A  “  Aj») 

•JaU^a)  »?r  l~nA 


Let  Ar  -  0  be  defined  as  the  position  at  which  »?r  *  0.5,  then: 

»?  a 

log  Tifr  -  -m* 

j-  »JA 


-A-  .  e"^A  ,  A - 1  -  e^A 

1“»»A  »»A 


•7A  “  i  + 

As  a  -*  +«,  ija  -♦  0,  i.e.,  the  original  flow  is  fully  depleted  by 
collisions  with  the  slow  atoms,  while  for  a  -*•  i,A  l,  as 
required  by  the  initial  conditions. 

with  the  above  formula  for  «;A'  the  normalized  slow  flow 
density  is 

*?s  -  n8/n0 

-  M(1  -  »7 a) 


+ 1)  [i  - 


iii±ii4A 

1  +  e^A 


wwv.v  v 
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As  A  -»  «,  »»A  •*  (m  +  1)  i  which  is  constant,  while  for  \  -  -®, 

0,  consistent  with  the  initial  conditions.  The  total 
density  at  any  position  is  then: 

IT  “  »»A  +  »»s 

.  I  ±  Qf+Dg"- 

1  +  e^A 

which  has  a  maximum  value  (p  +  1)  as  a  -  ®. 

Conservation  of  momentum  requires  the  development  of  a 
quantity  like  the  pressure  so  that: 

p  +  pu2  *  constant, 

where  p  -  mAn  is  the  mass  density.  The  pressure  change 
across  the  shock-layer  is  thus: 

Ap  -  ~A(pU2) 

n0UAz  "I 
(m+1)^ 

-  mAn0uA2  (^) 

Specification  of  the  pressure  change  would  determine  the  value  of 
p  for  a  given  n0uA2.  If  p  is  specified,  however,  conservation  of 
momentum  merely  demands  the  preceding  value  of  Ap.  Similarly, 
conservation  of  energy  will  determine  an  enthalpy  change  in  terms 
of  the  continuity  equation  and  specified  value  of  p  .  The 
conservation  principles  of  physics  thus  define  changes  in 
quantities,  such  as  pressure  and  temperature,  that  are  not 
necessarily  well-defined  in  a  non-equilibrium  flow  situation. 

The  transition-flow  shock-layer  as  formulated  resembles  a 
standard  shock  in  that  the  component  of  flow  velocity  normal  to 
the  layer  is  reduced  while  mass  flows  continuously  through  the 
layer.  The  characteristic  thickness  of  the  shock-layer  may  be 
estimated  by  the  distance  for  the  original  atomic  flow  density  to 
decrease  from  0.9  to  0.1  of  its  freestream  value: 


I. 


mA  LnouA 


2  _ 


(/i  +  1) 
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AA 


1 


^  (St? 


SLS\ 

o.i J 


i  log  9  -  ^ 


As  a  quick  example ,  the  density  ratio  across  a  hypersonic  shock 
in  an  ideal  (7  ■  5/3)  gas  is  (7+1/7-I)  ■  4,  so  n  »  3  and  the 
shock-layer  thickness  is  about  1.5  mean  free  paths. 

In  the  manner  of  an  oblique  shock,  it  is  possible  to  apply 
the  transition-flow  shock-layer  model  to  flow  onto  a  wedge.  The 
velocity  component  parallel  to  the  layer  is  preserved,  while  the 
component  normal  to  the  layer  is  reduced,  thereby  providing  an 
angular  displacement  of  the  flow  velocity  vector. 

us 

tan(0-6 )  -  — 

«A 


where  ut  is  the  flow  component  tangential  to  the  layer.  The 
relative  flow  speeds  normal  to  the  layer,  before  and  after,  are 
thus: 

uA  tan  p 

us  -  (m  +  1)  “  tan(0-«) 

With  the  distribution  of  flow  particle  densities  along  the 
direction  normal  to  the  shock-layer,  and  the  angular  arrangement 
of  the  wedge  flow  (Figure  32) ,  it  is  possible  to  consider  the 
distribution  of  atoms  ablated  from  the  wedge  surface.  For  the 
case  of  weak  ablation,  ablated  atoms  can  be  treated  as  impurities 
with  a  flow  density  nj  from  the  surface  that  is  diminished  by 
collisions  with  the  main  flow  atoms  (both  nA  and  ns) : 

dnj 

uj  -  -nj[nA  Qi  uIA  +  ns  Qj  uIs] 


where  ui  is  the  impurity  flow  speed  along  the  y-direction 
(outward  normal  from  the  surface)  and  Qi  is  the  cross-section  for 
collision  with  the  main  flow  atoms.  The  encounter  speeds  are 
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defined  by  the  geometry  of  the  particle  velocities: 

U1A2  "  uo2  +  UI2  “  2ucui  cos (w/2+6 ) 

and 

uifl2  ■  u82  +  u02  cob20  +  ui2  , 

where  uja  and  uj8  are  the  impurity  encounter  speeds  with  unslowed 
and  slowed  atoms,  respectively,  and  uQ  is  the  freestream  flow 
speed.  A  position  y  along  the  normal  to  the  surface  corresponds 
to  a  distance  along  the  normal  from  the  shock-layer  of: 

x  ■  i  sin(^-S)  -  y  cos (p—s ) 

where  i  is  distance  along  the  surface  from  the  wedge  vertex 
to  the  ablation  source.  In  terms  of  i  ■  n0Qx,  the  equation  for 
the  relative  impurity  density  »?x  “  ni/n0  is  then: 

^  dA~  "  cos(J-fi)  [q~)  +  *s(A>Ms] 

where  MA  -  uai/uj  and  Ms  -  uis/uI* 

With  *A  -  l/(l+e"A)  and  r,e  -  (M+l)mPx/  (l+e^A) , 
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fQl'l  1  l(A)  Ma  fl+e"As'| 

lojJ  co8^-5)  1o*  •  ma(a*“a)  +  —  Ml+sH 


,  fl+e^AS^| 


and  the  impurity  density  relative  to  its  surface  value  rjj  (is)  is! 

Ma-(m+1)Ms  Qi/Q 

»»l(A) 


»»l(As) 


ri+e*i*s'i 

ll+e^A  J 


e“MA(AS“* ) 


COS  (0-5  ) 


This  may  be  rearranged  for  better  inspection  as: 


MA-(p+l)M8 

Ql/Q 

»7l(A) 

ri+e  *  -ru+i)MB(Afl-A) 

COB (0-6 ) 

»?l(As) 

L+«V  J  e 

.  4 

► 

t 


which  indicates  a  modified-form  of  exponential  decay  away 
from  the  surface. 

For  A  -  0,  and  e">As  «  1, 


il(0) 

m(As) 


[Ma-C^+DMsIQ!/© 


>iCOS  (0-6) 


-Hm+1)  MsQ!/Q] 

COS (0-5 ) 

e 


Since  As  -  n0Qi  sin(0-«)/  the  impurity  level  at  the  shock-layer 
should  decrease  exponentially  with  slant  height,  i.  This  result 
would  be  true  only  if  e”^As  was  small  for  all  values  of  i.  Near 
the  vertex,  however,  impurities  can  reach  the  relatively  low 
density  regions  of  the  main  flow  more  readily.  These  impurities 
can  then  be  convected  along  lines  between  the  shock-layer  and  the 
surface,  arriving  at  distances  from  the  surface  that  could  not  be 
attained  by  normal  diffusion  through  the  higher  density  region 
near  the  surface. 

It  is  possible  to  estimate  conditions  near  the  leading  edge 
of  the  flow  deflection  by  adapting  the  shock-layer  analysis.  If 
a  shock-layer  at  angle  0  to  the  freestream  is  extended  toward 


the  vertex  of  the  wedge,  the  distance  between  the  shock-layer  and 
the  wedge  surface  will  become  insufficient  to  allow  the  flow  to 
turn  through  the  required  angle  s .  Flow  particles  will  therefore 
impinge  on  the  wedge  surface  and  must  be  re-emitted.  At  any 
point  between  the  shock-layer  and  the  surface  the  net  flow  will 
comprise  populations  of  particles  moving  in  different  directions. 
To  estimate  conditions  near  the  leading  edges,  these  populations 
must  be  modeled. 

In  the  previous  analysis,  all  particles  not  traveling  with 
the  velocity  of  the  unscattered  particles  (uA  normal  to  the 
shock-layer)  were  assigned  a  (normal)  speed  u8  and  a  density  ns. 

A  return  flux  of  particles  from  the  surface  has  the  effect  of 
reducing  the  speed  u8.  It  may  be  reasonable,  therefore,  to 
distribute  the  total  particle  flux  n8u8  between  two  populations: 

n8u8  «  n8iu8i  -  nS2Us2 

where  the  subscripts  '1'  and  '2'  refer  to  particles  flowing 
to  and  from  the  surface,  respectively.  The  actual  distribution 
of  particle  speeds  and  densities  may  be  expected  to  vary  within 
the  broad  constraint  of  total  flux,  n8u8,  as  a  function  of 
position  between  the  shock-layer  and  surface.  It  is  only  at  the 
surface  itself  that  an  additional  constraint  can  be  applied, 
namely  that  of  zero  net  flux  to  the  surface.  Without  invoking 
additional  conditions  associated  with  the  energetics  of  the 
particle-surface  collision,  the  only  characteristic  speed  in 
the  problem  is  the  initial  normal  component  uA.  It  is  useful, 
therefore,  to  require  one  population  to  have  the  speed  uA  and  to 
solve  for  the  necessary  speed  of  the  other  population.  Thus,  let 
usl  “  UA  and  us2  *  uw  such  that: 


n8u8  -  njuA  -  n2uw 


From  the  geometry  sketched  in  Figure  33,  uA  =  uc  sin/9,  where  uc 
is  the  freestream  flow  speed.  At  the  surface,  the  particle  flux 


[nA(A8)  +  ni(As)  ]uA  sins 


n2uW 
COS (0-6 ) 


ft 


Proa  the  previous  analysis,  we  have  the  condition  on  the 
nornal  particle  flux  through  the  shock-layer: 

no«A  “  («A  +  «l)UA  “  n2uw 

Furthermore,  the  average  speed  of  the  slowed  particle  flow  is 
given  by: 


«A 

u„ 


tan/9 


tan (p-6 ) 


so 


u  .  u.  fcanLg-iJ 

u®  UA  tan/9 


Combination  and  re-arrangement  of  equations  provides  the  speed  of 
particles  returning  from  the  surface  (in  the  direction  normal 
':o  the  shock-laver) : 


where  17 A  is  n^/no  for  x  -  \B. 

The  relative  density  of  this  flow  population  is: 

n2  uA  sini  cos (p-6) 

nQ  "  uw  [sin/9  -  sim  cos  (/9 -6)] 

(The  relative  density  of  the  other  population  is: 
ni  sin/9 

nQ  ”  [sim  -  sim  cos(^-s)  ]  “  ' 

with  fjA  again  evaluated  at  the  surface.) 

In  the  earlier  analysis,  when  the  distance  of  the  shock- 
layer  from  the  surface  is  sufficient  to  establish  a  single  slow- 
particle  population  traveling  parallel  to  the  surface,  atoms  from 
the  freestream  primarily  collide  with  these  slow  particles  and 
are  largely  prevented  from  reaching  the  surface.  The  preceding 
analysis  allows  a  transition  to  the  regime  near  the  leading  edge 
in  which  freestream  particles  can  bombard  the  surface  and 
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collisions  occur  bstwssn  re-emitted  particles  and  the  freestream. 
The  necessary  speed  of  particles  from  the  wall  can  be  calculated 
within  the  structure  of  these  analyses, 


uw 


Ue  "  cos (0-6) 


It  should  be  noted  that  the  angle  between  the  shock-layer 
and  the  freestream  flow  need  not  be  constant  as  the  wedge  vertex 
is  approached.  The  value  of  0  measured  at  distances  well  down¬ 
stream  of  the  leading  edge,  where  it  is  convenient  experimen¬ 
tally,  is  not  necessarily  the  value  needed  in  the  equation  for 
U£.  Indeed,  the  speed  of  re-emitted  particles  is  actually 
subject  to  an  additional  constraint  due  to  the  boundary  condition 
on  surface  heating.  For  example,  an  adiabatic  wall  (zero  heat 
transport  into  the  material  in  steady  state)  requires  the  same 
relationship  of  U£  to  the  speed  of  the  impinging  flow  as  in  free 
molecular  flow,  independent  of  the  flow  density: 


3-|( 


Si) 


E1LS1 


4* 


[i  ♦  (£?)/•*]} 


where  7  is  the  specific  heat  ratio  for  the  freestream  flow, 
S2  is  (7/2)  times  the  square  of  the  freestream  Mach  number  and 
the  modified  recovery  factor  for  an  adiabatic  wall  in  free 
molecular  flow  is: 


r’(S) 


S2 


( 


2S^  +  1  - 


1  +  7*(S  sins)  El  +  erf (S  sinfi)]e(S  sin5)  J 


(From  Schaff  and  Chambr6,  Flow  of  Rarefied  Gases.  Princeton 
University  Press,  1961,  p.  15.) 

Thus,  ue/u0  is  determined  by  surface  boundary  conditions  (allow¬ 
ing  ablating  or  conducting  walls  to  provide  different  flow 
results  than  insulating  walls)  so  the  local  shock-layer  angle  p 
can  be  calculated  and  used  in  the  other  formulas  in  a  self- 
consistent  fashion.  In  order  to  accomplish  this  accurately, 
the  shock  shape  must  be  constructed  relative  to  the  wedge,  so 


.that  i?a(as)  can  ba  obtained.  This  may  be  done  by  fixing  U£/u0 
and  6  and  then  varying  0  to  generate  ^a(ab)  “  1/(1  +  e^As) .  The 
value  of  xu(0)  provides  the  distance  n0Qxs  normal  to  the  shock- 
layer  (at  local  angle  0  to  the  freestream)  between  the  surface 
and  the  shock-layer.  (Note  that  it  is  not  necessary  for  x  -  0, 
corresponding  to  ij a  ■  0.5,  to  lie  beyond  the  wedge  surface,  since 
the  shock-layer  actually  extends  with  i»a  <  1  to  x  -  -  ®.) 

If  it  is  postulated  that  energy  accommodation  between  the 
flow  particles  and  the  surface  is  perfect,  then  the  thermal 
processes  re-emitting  flow  particles  also  drive  the  ablatants  and 
we  may  expect  that  the  ablatant  speed  is  related  to  flow  particle 
speed  in  terms  of  the  relative  atomic  masses  mx/mAt 


This  value  of  impurity  speed  can  be  used  in  the  analysis  of 
impurity  diffusion  from  the  surface,  indicating  that  lower  mass 
ablatant  atoms  can  penetrate  further  from  the  surface  since  the 
values  of  Ms  and  Ma  are  reduced;  the  value  of  Ms  is  also  less 
near  the  leading  edge  because  the  effective  value  of  u8  is 
reduced  by  the  re-emitted  flow.  (Lower  values  of  us,  however, 
increase  (n  +  1) ,  so  the  principal  gain  in  penetration  is  merely 
the  increased  proximity  of  lower  density  flow  closer  to  the 
vertex  of  the  wedge,  as  the  absolute  distance  between  shock-layer 
and  surface  decreases  relative  to  the  freestream  mean  free  path, 

l/n0Qi.) 

As  impurities  originating  near  the  vertex  convect  downstream 
along  the  lower  density  portion  of  the  shock-layer,  their  density 
will  diminish  as  diffusion  against  the  free  stream  and  back 
toward  the  surface  occurs.  Such  diffusion,  however,  is  limited 
by  the  short  mean  free  path  relative  to  the  increasing  gap 
between  the  shock-layer  and  the  surface.  Thus,  while  the 
relatively  high  density  main  flow  near  the  surface  can  restrict 
the  migration  of  ablated  material  from  most  of  the  surface,  it 
does  not  protect  the  surface  near  the  leading  edge.  (In  fact, 
ablated  material  off  the  leading  surfaces  of  flow  deflection  is 


prevented  from  participating  in  any  ameliorating  recollection/ 
equilibrium  with  the  later  surface) .  This  suggests  that  the 
greatest  points  of  erosion  in  high  speed  flow  channels  will  be  at 
the  initial  portions  of  flow  deflection  surfaces,  where  the 
distance  between  the  "shock"  responsible  for  deflection  and  the 
surface  is  less  than  a  few  mean  free  paths.  (Viscous  boundary 
layer  growth  can  result  in  deflection  of  the  free  stream  flow  so 
that  the  surface  within  a  few  mean  free  paths  of  the  leading  edge 
of  a  plate  at  zero  angle  of  attack  will  also  be  relatively 
unprotected.) 


The  theoretical  formulation  of  the  preceding  section  can  be 
used  to  examine  experimental  results  from  the  spectroscopic 
studies.  In  utilizing  the  formulas,  it  is  necessary  to  evaluate 
the  appropriate  shock-layer  angle,  p,  for  the  location  of 
interest,  especially  near  the  leading  edge.  For  example,  by 
increasing  p  at  fixed  deflection  angle  s ,  it  is  found  that  the 
condition  p  -  90*  corresponds  to  r/A  m  1»  and  thus  represents  the 
leading  edge  of  the  flow  deflection.  At  6  ■  22*,  however,  va  K  0 
for  values  of  p  <  60*,  a  situation  that  indicates  sufficient 
distance  between  the  shock-layer  and  the  surface  to  turn  all  of 
the  freestream  atoms  through  the  necessary  deflection  6 .  Lower 
values  of  p  represent  a  shock-layer  separated  from  the  surface  by 
a  collisional  flow,  with  the  value  of  p  provided  by  the  usual 
oblique  shock  relations.  The  formulation  as  constructed  thus 
provides  a  transition  from  free  molecular  to  continuum  flow  as 
the  shock-layer  curves  from  normal  angle  of  attack  at  the  leading 
edge  to  a  constant  shock  angle  a  few  mean  free  paths  further 
downstream. 

In  the  immediate  vicinity  of  the  leading  edge,  it  is 
possible  to  estimate  the  diffusion  of  surface  ablatants  into  the 
flow.  For  example,  if  we  consider  the  experimental  situation 
of  a  partially  ionized  argon  flow  (y  =  1.4)  onto  a  wedge  with 
S  «  22*,  the  surface  re-emission  speed  relative  to  the  freestream 
speed  is  ue/uq  *  0.26.  At  a  shock-layer  angle  of  p  -  85*,  the 
calculated  density  of  freestream  atoms  impinging  on  the  surface 
relative  to  the  freestream  value  is  ija  =  na/n0  *  0.56.  This 
relative  density  corresponds  to  a  normalized  position  in  the 
shock-layer  As  -  -0.079.  That  is,  the  wedge  is  actually 
penetrating  slightly  more  than  half  way  into  the  transition 
flow  shock-layer. 

The  diffusion  of  surface  atoms  occurs  at  relative  speeds 
normal  to  the  surface  given  in  terms  of  the  mass  ratio  of 
ablatants  compared  to  the  argon  atoms  of  the  re-emitted  flow. 
Table  II  provides  relative  speed  values  for  carbon,  fluorine,  and 


hydrogen,  along  with  the  densities  of  these  impurities  relative 
to  their  surface  values  for  a  normalized  distance  a  -  -1  (in  the 
direction  normal  to  the  shock- layer) .  A  principal  factor  in 
determining  the  relative  densities  is  the  relative  cross-section 
Ql/Q,  so  different  values  are  indicated  in  Table  II.  In  partic¬ 
ular,  if  Qj/Q  -  1  is  used  for  both  carbon  and  fluorine,  the  lower 
atomic  mass  of  carbon  allows  it  to  penetrate  further  than  the 
fluorine.  This  results  in  more  than  a  fifty  percent  higher  value 
for  carbon  yg  fluorine  compared  to  expected  values  based  on  a 
surface  chemistry  for  Teflon  (CF2)  that  would  favor  a  higher 
fluorine  density  in  the  flow.  If  the  cross-section  for  carbon- 
on-argon  is  smaller  than  fluorine-on-argon,  the  relative  density 
of  carbon  is  enhanced  further,  as  indicated  by  a  calculation 
using  Qj/Q  *  0.5  for  which  the  relative  carbon  density  is  22 
times  that  of  fluorine.  (Experimentally,  carbonization  of  the 
surface  with  repeated  test  operation  may  cause  spurious  results. 
The  leading  edge,  however,  tended  to  remain  free  of  such  carbon.) 
As  a  separate  comparison,  hydrogen  (for  a  polyethylene  surface) 
has  a  higher  relative  density  than  carbon  due  to  both  its  lower 
atomic  mass  and  a  possibly  smaller  cross-section  (Qj/Q  =  0.67) 
relative  to  the  argon-argon  collisions  determining  the  main  flow 
field. 

For  a  freestream  density  of  n©  ■  1022  m"3  and  a  cross- 
section  Q  -  25  x  10“20  m2,  the  characteristic  mean  free  path  that 
scales  the  flow  near  the  leading  edge  is  i  “  0.4  mm.  Experi¬ 
mental  data  indicate  that  carbon  and  fluorine  are  present  with 
fairly  uniform  distributions  (vs  strong  exponential  decay  from 
the  surface)  for  distances  of  0.2  to  0.7  mm  from  the  surface  for 
fluorine  and  0.2  to  0.9  mm  for  carbon  at  a  lateral  distance  of 
2.5  mm  from  the  wedge  vertex.  At  larger  lateral  distances,  the 
concentrations  of  carbon  and  fluorine  are  actually  greater  at 
several  tenths  of  a  millimeter  into  the  flow  than  at  the  surface, 
with  a  thickness  (normal  to  the  surface)  of  a  few  tenths  of  a 
millimeter.  Since  the  mean  free  paths  for  penetration  of  the 
main  flow  at  distances  well  back  from  the  leading  edge  are  much 
less  than  the  distance  to  these  concentrations,  it  appears  that 
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TABLE  II 


6  -  22*,  fi  -  85*, 


-1.0 


Material 

ui/u0 

Relative 

Speed 

Ql/Q 

Relative 

Cross-Section 

ni/nxw 

Relative 

Density 

Carbon 

0.467 

1.0 

1.9  x  10“2 

0.5 

1.4  x  10"1 

Fluorine 

0.371 

1.0 

6.1  x  10“3 

Hydrogen 

1.62 

0.67 

2.1  X  10”1 

this  material  has  convected  from  the  more  diffusive  region  within 
a  mean  free  path  of  the  wedge  vertex. 

While  considerable  uncertainty  exists  concerning  the  proper 
values  to  use  for  quantities  such  as  collision  cross-sections  for 
the  chemically-active,  non-equilibrium  plasma  flow,  the  present 
theoretical  formulation  could  be  applied  to  the  modeling  of 
plasma/surface  interactions  in  actual  devices.  For  such  applica¬ 
tion,  it  would  be  useful  to  benchmark  the  collisional  processes 
in  a  separate  apparatus  with  appropriate  diagnostic  access.  In 
this  way,  the  relevant  data  base  for  the  materials  and  particle 
energies  of  interest  for  a  particular  project  could  be  obtained. 
Even  without  an  accurate  absolute  data  base, scaling  between 
device  operation  and  the  idealized  sample  tester  could  be 
performed  in  the  context  of  the  theoretical  analysis,  so  that 
estimated  component  lifetimes  in  high  specific  power  devices 
could  be  made. 
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